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O Abstract 

(N ■ 

In this paper we derive the one-dimensional bending-torsion equilibrium model modeling the junction 
of straight rods. The starting point is a three-dimensional nonlinear elasticity equilibrium problem 
written as a minimization problem for a union of thin rod-like bodies. By taking the limit as the 
thickness of the 3D rods goes to zero, and by using ideas from the theory of T-convergence, we obtain 
that the resulting model consists of the union of the usual one-dimensional nonlinear bending-torsion 
rod models which satisfy the following transmission conditions at the junction point: continuity of 
displacement and rotation of the cross-sections and balance of contact forces and contact couples. 

< 

en ■ 1 Introduction 

> | 

■ In many real-life situations, such as, for example, in certain types of bridges or building structures, two (or 
I, , several) elastic rods are connected at one point. Such points where several rods meet are called junctions. 

Such multiple rods systems may be as small as two rods joining in a non-smooth way, or as complex as 
' several hundreds of interconnected rods forming a massive network. In either case, the basic principles 

■ of analysis are the same (although the complexity of the computation depends on the complexity of the 
system). Therefore, in the present paper, we limit our study to the case of one junction point. 

| In this paper we consider the equilibrium problem of a three-dimensional elastic body which consists 
of n straight thin rod-like bodies connected in a single point. Since the rods are thin, the behavior of each 
^) , rod should be well approximated by the one-dimensional rod model. In order to close the model one needs 
,— 1 \ conditions at the junction point. These conditions can be seen as transmission conditions as well. Since we 
T 7"! ■ are interested in the bending-torsion behavior of rods, such a rod is expected to be governed by the fourth 
^ . order equation, see 0j. Since this equation can be written as a first order system in terms of the contact 
^ \ force, the contact couple, the rotation of the cross-section and the deformation (displacement), we expect 
^ ■ the following four junction conditions (based on the continuity of the deformation and equilibrium laws): 

1) the sum of all contact forces at the junction is zero, 

2) the sum of all contact couples at the junction is zero, 

3) continuity of the rotation of the cross-section (the angles at the junction point are preserved), 

4) continuity of the displacement (deformation/position) at the junction point. 

These conditions follow physical intuition and are already used in modelling networks of elastic rods (see 
e.g. [13] and [31]; in the case of strings see [TT]). 

In [14] these junction conditions have been mathematically justified for the case when the starting 
configuration is that of three-dimensional linearized elasticity. We justify these junction conditions starting 
from a three-dimensional nonlinearly hyperelastic material in the formulation as the energy minimization 
problem. Since the rods are thin, we recognize the small parameter h describing their thickness. The 
mechanical response of rods strongly depends on the relative magnitude of the applied load with respect 
to the rod thickness h. In [23] a bending-torsion model of a single nonlinearly elastic inextensible rod was 
derived by the theory of T-convergence and the geometric rigidity theorem from [12]. In order to obtain 
the bending-torsion model the main assumptions is that the energy of the rod is of the order /i 4 . For other 
models see [U [26]. In the present paper we would like to obtain junction conditions at the junction of 



rods for the case when the total energy functional is of order /i 4 . However, unlike in the case of a single 
rod studied in [12], in the case when 2 or more rods meet at a junction, we cannot rescale our problem in 
such a way that the entire problem is defined on a canonical domain independent of h, at least in a simple 
way. To deal with the complications related to the geometry at the junction, we assume that the junction 
region of the rods forms a domain which scales with h (say a sphere, at which all the rods are connected). 
Then, as h — > 0, the junction region converges to a point. This leads to a problem with no obvious simple 
canonical domain, and so the results from [24] cannot be applied directly to this problem. To get around 
this difficulty we adapt the ideas from |24] to this new scenario and express the asymptotic behavior of 
minimizers in norms depending on the thickness h. 

Following |24] , we first prove a compactness result (Theorem Theorem I3.1|) for the sequence of energy 
minimizers y^ 1 ' deriving the asymptotic behavior of Vj/^*'. Moreover we prove that the rotations of the 
cross-sections need to be continuous at the junction point in the limit as h — > 0. Since we are considering 
a pure traction problem for rods joining at a point, we still need to control the displacement of the entire 
structure. Under the assumption that the translation of the whole structure is controlled at the end of 
one rod, in Lemma 14.11 and Corollary 14.21 we derive the asymptotic behavior of the minimizing sequence 
yV 1 ' and we obtain that in the limit, the displacement (deformation) of the rods at the junction point is 
continuous. Finally, in Theorem 15.11 we derive the model for the junction of rods. 

Junction of elastic rods has been studied by several authors. However, most results are restricted to 
linearized elasticity. The first study of the junction of two rods is given by Le Dret in [19J, see also [21J and 
[30] . For systems of rods see also [27] and [28] and references therein. The junction of two plates is studied 
in [15] and [201 E2], while [29] deals with the junction of beams and plates. The case of the junction of a 
three dimensional domain and a two dimensional one is explored in [10], see also [9] and references therein. 
For the asymptotic analysis of the junction between three-dimensional structures and one-dimensional one 
see [18] and [5]. See also 0[23] for the asymptotic analysis of the problem of junctions of thin pipes filled 
with a fluid using asymptotic expansion method. 

Two efforts in the study of junction problems within nonlinear elasticity are made in [6] and [16] using 
asymptotic expansion method. In [6] the model of plate inserted in a three-dimensional elastic body is 
derived, while in [16J a model of junction of rod and plate is derived. 



2 Setting up the problem 

The domain of the junction of rods we define as a union of cylinders and the "junction" part. Let n £ N 
denote the number of rods meeting in junction and let h > 0. Let each rod be of length Lj with the cross- 
section hSi, where Si C M 2 (open, bounded, connected). Let the junction part is of the form T h = hT, for 
T C M 3 open, bounded, connected set. Let Q i £ SO(3), i = 1, . . . ,n. The vector tj = Qjei denotes the 
tangential direction of the i-th rod. Then the domain of the junction of rods is given as 

n 

n h = T h U (J Ct C } : = Qi((h, U) x hSi). 

1=1 

We assume that the domain £l h is open, bounded, connected and with the Lipschitz boundary. We also 
assume, as in [21], for each i that 

/ X2XsdX2dX3 = / X2dX2dX3 = / x^dx2dxs = 0. 

Every function y G W l,2 ((a, b); M. 3 ) we naturally interpret as an element of H /1,2 ((a,6) x M 2 ;R 3 ). We 
also define the mapping : (a,b) x Si — > (a,b) x hSi by P^(ii,i2,X3) = (x±, hx2, hx^) and use it to 
change between thin and thick domain. 

The starting point of our analysis is the equilibrium problem of the junction of rods, i.e. elastic body 
i} h . The internal energy of the junction of rods is given by 

E {h \y) = f W(Vy(x))dx, 

for a deformation y G W 1 ' 2 (fl A , R 3 ), where W : M 3x3 — > [0, +oo] is an internal energy density function. 
For W, as in [24], is supposed to satisfy 



• We C°(M 3x3 ), W is of class C 2 in a neighborhood of S0(3); 

• W is frame-indifferent, i.e., W(F) = W(RF) for every F G M 3x3 and R G SO(3); 

• W(F) > C w dist 2 (F, S0(3)), W(F) = if F G SO(3). 

We are looking for the one-dimensional b ending-torsion model of junction of rods. Thus, motivated by 
|24| . we will assume that the energy E^ 1 behaves as h A . Then we analyze the behavior of £^'(j/)//i 4 and 
derive the one-dimensional model. This is in [24J obtained by T-convergence, but in the junction problem 
there is no obvious and simple canonical domain, a domain that is independent of the thickness h. Still, 
using the ideas and techniques of T-convergence we are able to give the asymptotics (in the form (|4.17|) ) 
of the infimizing sequence of the total energy functional and the total energy functional itself. 
We shall need the following theorem which can be found in |12j . 

Theorem 2.1 (on geometric rigidity) Let U C R m be a bounded Lipschitz domain, m > 2. Then there 
exists a constant C{U) with the following property: for every v G W 1,2 (U; M m ) there is associated rotation 
R G SO(m) such that: 

\\Vv-K\\ L 2 (u) < C(U)\\di S t(Vv,SO(m)\\ LHu) . (2.1) 

We will apply this theorem in the next section on subdomains of Q h which are of size h in each direction. 
This is possible since the constant C(U) in the estimate is independent on the translation and dilatation of 
U. Let us consider the domain hll, for h > 0. Take v G W 1,2 (hU; R m ). Then the function 
belongs to W 1,2 (U; M. m ) and satisfies the estimate 

\\Vv^-R\\ L 2 {u) <C(U)\\dist(Vv^,SO(m)\\ L 2 {u) . 

Since \7v^ = X7v(hx) after the change of variables in the norms we obtain that the estimate (|2.1|) holds 
for v with the same constant C(U). See also [12J. 

Throughout the paper we use the following function space 

iy 1 ' p (0;SO(3)) = {R G ^ 1 ' p (fi;M 3x3 |R(x) G SO(3) for a.e. x G O}. 

Moreover, by || • || (without subscript) we denote the Frobenius matrix norm. 



3 Compactness 

In this section, following [23] we prove the compactness result (Theorem 13. ip . Namely, for that satisfy 
()3.ip (this will be shown for infimizers j/W of the energy of order /i 4 ) we obtain asimptotics of X7y( h h 
Moreover it turns out that rotations of the cross-sections in the limit, when h tends to 0, need to be 
continuous in the junction point. 

Theorem 3.1 Let (yW) C W l > 2 (tt h ;R 3 ) be such that 

limsup-^ / dist 2 (Vy W ,SO(3))cfe < +oo. (3.1) 

Then there exists a subsequence (not relabeled) and Rj G W 1,2 ((0, Li), SO(3)), i = l,...,n such that 
Ri(0) = R2(0) = • • • = R n ,(0) in the sense of traces and 

1 n f — 
lim T2E/ \\Vy {h) (x)-R l (x-t i )\\ 2 dx = 0. (3.2) 



Proof. We follow proof of Theorem 2.2. in 

Now we cover £l h with subdomains of size h in each direction and apply Theorem 12.11 on each of them. 
For every h > and i = 1, . . . , n let G N be such that h < L{/k^ < 2h and let 

r ak y.= (a,a + ^\, ae[0, Ll )n^N. 



We apply Theorem 12.11 to domains Qj((a,a + 2h) x hSi) (when a = Li — rfc we take (L% — 2h,Li)) and 
T h UUl l =1 Qi({h,2h) x hSi). Note that Q;(P fc x hSi) C Qi((a,a + 2/i) x hSi). Then there exist a constant 

a,K i 

C (independent of % (as there is finite number of domains) and h (by a note after Theorem 12. ip ) and 
a piecewise constant map : U™ =1 Qj([0, Li] x {0} x {0}) — > SO(3), constant on each [a, a + -^) for 

a G [0, Lj) n and on T h U U" =1 Qj([/i, r£) x /iS^), such that for every i 6 {1, . . . , n} we have: for every 

ae[o,L,)nf|N 

/ ||Vy w -R (/i) || 2 dx < C [ dist 2 (Vy (/i) ,SO(3))dx 

•'QiU i h xhSi) JQi((a,a+2h)xhSi) 

a,k. 
i 

and 

/ ||Vy (h) -R (/l) || 2 cfx < C / dist 2 (Vy (/l) ,SO(3))(fx. 

i 

By summing all these estimates, since only neighboring subdomains overlap, we obtain the inequality 

A / l|Vy W -K {h) \\ 2 dx < ^ [ dist 2 (Vy (/l) ,SO(3))dx < dh 2 , (3.3) 
nr Jnh h z J nh 

where the last inequality holds for h small enough by (|3.ip . 

In the sequel we show that on a subsequence R^ converges to a W 1 ' 2 function. In order to do that 
we first estimate the difference of R^ 1 ) on neighboring subdomains. 

Let now a% G (0, Li — 4/i]n If N, foj = Ui+yt ■ Now we apply Theorerr l2.il on the set Qi((cij, aj + 4/i) x /iS^) 

i i 

we obtain that there exist R £ SO (3) such that 

f \\Vy {h) -R\\ 2 dx < C 2 [ dist 2 (VyW,SO(3))da;. 

jQidat ,a,i+4h) xhSi) J Qi{(ai,ai+4h)xhSi) 

Then using the facts that I* k h^\. k h are contained in (dj, o« + 4/i) x /i5j we have for every i: 
^||RW(M i )-RW(M i )|| 2 
< 2§ (||RW (aiti) - R|| 2 + ||R - R (h) (6i*i)|| 2 



< 2 / ||R (h) M 4 ) - R|| 2 + 2 [ ||R - R^(6A)|| 2 

jQiCJ* uh xhS t ) JQ t (P uh xhSi) 

<4 / ||RW( ait .)_VyW|| 2 + ||Vy (h) -R|| 2 



+ 1 / || R - V-y""|| J + \\Vy {h) - R W (M; 

.CO _WI|2 



Qi(J« .xMi 



<4/ ||R( h )(aiti) - Vy^|| 2 + 4 /" ||Vy^ - R|| 

+4/ ||VyW-RW(6A)f- 



All terms on the right hand side of the estimate can be estimated by Theorem 12. 1\ so we obtain 

^||R(»)(aA) -R^Mi)!! 2 < § /" dist 2 (VyW,SO(3))dx, (3.4) 

fc i " JQ i ((o i ,o<+4fe)x/iS i ) 

and similarly, as 7^. are contained in T h U Qj((/i,4/i) x hSi), we obtain 

77t> i 



k 

^IIRW(O) -RW(§t 4 )|| 2 < § / dist 2 (VyW,SO(3))dx. (3.5) 



Thus we have (since R^ is piecewise constant) for every < £ < -rfc and every i and for every a G 
(0, U) n f£N s.t (a, a + Ah) C (0, Li): 

/" \\K^((x 1 +0t i )-^ h \x 1 t i )\\ 2 d Xl < S /" dist 2 (Vy^ SO(3))dx, (3.6) 

t " yQi((a,a+4/i)x/i.5i) 



since xi + £ and x\ belong to the same or neighboring subdomains and we can apply estimate (|3.4p . In 
the same way we can show that for every i and a s.t. (a — 2h, a + IK) C (0, Li) and every — Lj/fcf < £ < 0, 

[ WR^dxi+Ot^-n^Hxi^fdxi < § /" dist 2 (Vy^,SO(3))dx. (3.7) 

J/ 1 « jQ l ((a-2h,a+2h)xhS i ) 

7 z 

Let us now look at cylinders C\ and C 2 . By summing estimates (13.4p . (|3~5j) . (13.6p . (13.7p we have that 
for every open interval I' compactly contained in (— Li,L 2 ) and £ G M. which satisfies for all i, |£| < 
dist(/',{-L l5 L 2 }),|el < f 

/ HR^^i+a-R^^i)!! 2 ^! < £ / dist 2 (VyW,SO(3))dx, (3.8) 
Jl> h z J nh 

where R^ } : (-L 1 ,L 2 ) SO(3) is defined by 

B (% T x _ / R (h) ("Mi), ifa;iG(-Li,0] , , 

RW(xi t 2)) if Xie(0 ,L 2 ) • (3 - 9) 

By iterative application of (|3.8p and using the inequality (xi + • • • + x n ) 2 < n(x 2 + • • • + x 2 ) and the 
assumption (|3.ip for every open interval /' compactly contained in (— Li,L 2 ) and f 6l, which satisfies 
\i\ < dist(/',{-Li,L 2 }), we have 

J f ||RW(X! + - RW(x 1 )|| 2 dx 1 < C 4 P-Q + lY p J dist 2 (VyW, SO(3))dx < C 5 (|£| + /*) 2 , (3.10) 

Note here that the factor (^ + l) 2 is the upper estimate of the number of terms by which the left hand 
side of (|3.10p has to be estimated. Using the Frechet-Kolmogorov (see [2j Theorem 2.21, Theorem 2.22]) 

Q 1 )! 2 ) \ 

criterion, one can deduce from this that for any sequence hj — > there exist a subsequence (R m ' ) 
strongly converging in L 2 {— Li,L 2 ) to some R G L 2 (— Li,L 2 ) with R(xi) G SO(3) for a.e. x\ G (— Li,L 2 ). 
We define R~i : (0,Li) -> SO(3),R 2 : (0,L 2 ) -> SO(3) as 

Ri(xi) = R(-xi), if xi G (-Li,0), 
R 2 (xi) = R(xi), ifxiG(0,L 2 ). 

We shall prove that R G W 1 ' 2 ((-L 1 , L 2 );M 3x3 ). Using the estimate (pHUT) and letting h -> we obtain 
that for every /' compactly contained in (— Li,L 2 ) and every £ which satisfies |£| < dist(I', {— L\, L 2 }) 
there exists a constant C independent of /' and £ such that 

|R(xi +0 -R(xi) 112 



-dxi < C. (3.11) 



From standard theorems we obtain that R G W 1,2 ((— L\, L 2 ); M 3x3 ). This is equivalent to the fact that 



L_ _ _ _ 

Ri G TU 1 ' 2 ((0,Li);lR 3><3 ), R 2 G iy 1 ' 2 ((0, L 2 ); M 3x3 ) and Ri(0) = R 2 (0) in the sense of traces. In the same 

way one can take cylinders C\ and for i = 3, . . . , n (by choosing every time a subsequence R^ 1 '-'* of 

already chosen sequence R 1 >— ) we have the existence of R^. Moreover, the definition of Ri is not 
ambiguous and Ri(0) = R 2 (0) = • • • = R n (0). 

Ttitf l|Vy (/l) (x)-R i (x-t,)|| 2 dx<-| f \\VyW-RWfdx + ^ E f W^i* ■ U) - R {h) (x)\\ 2 dx 

Using the estimate ()3.3p and RW — >• Rj in L 2 (0, Lj) we obtain that 

1 n r — 
lim n V / ||VyW(x) - Ri(x • *i)|| 2 d& = 0. 



4 T-convergence 

In the Theorem 13.11 we obtained the asymptotics of X7y( h \ Still, as we are in the pure traction case, in 
order to obtain the asymptotics of y^ one needs to control the constant. Thus we additionally assume 
that the mean value at the end of the first rod behaves nicely. Then we obtain that in the limit in the 
junction point displacements from different rods have to be equal. 

Lemma 4.1 Let (hi) be a sequence that converges to and (y^ h ^) C W 1,2 (£l hj ; R 3 ) such that 

lim sup — 2 / ||Vy (h ^|| 2 dx<oo. (4.1) 
j->oo h- J < >'• , 

Let there exist y® £ W 1,2 ((0, Li); M 3 ) such that y®(0) = in the sense of traces and let us suppose that for 
every i, 

lim f l|V(y^ o Qi o _ ( ( y oy | q | )\\ 2 dx = 0. (4.2) 



Let us also suppose that there exists 



lim f yte) o Q x o pCy) dx := C Ll G 



(4.3) 



Then for Co := Cl 1 — y\(L{) we have 

lim o Q x o pfy) - C \\ L 2 ({h]}xSl) 



J->00 



lim ||y^ o Q„ o P^') - C || L 2 ({MxSn) = (4.4) 



J-5>00 



and 



lim £ o Q, o P(^) - || w i,2( (fti ,L 4 )x5 1 ) = 0, 



1=1 



Proof. By applying the Poincare inequality (see part b) of [U Theorem 6.1-8]) to the cylinders (0, 1) x 
we have that there exists a constant K\ such that for every i £ {1, . . . , n} and every y £ W 1,2 ((0, 1) x S 1 ,, 
one has 



y - -f ydx < Ki\\Vy\\L2((a t i)xSi^)- 

^ xSi L 2 ((0,l)xS i; R 3 ) 

By applying this estimate on functions of the form y(x) = y((Li — h)x% + h,X2,x$) we obtain that there 
is a constant K% = max{l, Lj — h}K\ such that for all i € {1, . . . , n} and for all h > (small enough) and 
all y £ W^ 2 ({h,L t ) x S^R 3 ) one has 



y — x ydx 

{Li}xSi 



< K 2 \\Vy\\ L 2^ h;Li)xSi . ]l 



(4.5) 



In a similar way we obtain 



y — x 

{h}xS l 



< K' 2 \\Vy\\ L 2 { ( hM) 



xSi 



(4.6) 



L 2 ((h,L l )xS l -l 



Moreover, by using the same rescaling of the domain, from continuity of traces we obtain that there is a 
constant K3 such that for all i,h and y £ W ' 2 ((h, Li) x S'jjM 3 ) one has 



\y\\L 2 ({h}xS l -M 3 ) + \\y\\L 2 ({Li}xSi;R 3 ) - ^WvWw^dhM) 



xSi 



(4.7) 



By applying the Poincare inequality (of the same form as before) to the domain T on functions given 
by y(x) = y(hx) we have that there exist a constant K4 such that for all i,h and y G W l,2 (T h ; R 3 ) one 
has 



y — X ydx 

Qi({h}xhSi) 



< hK4\\Vy\\ L 2r T h. 



L 2 (T h -J 



In the similar way as before we conclude that there exists a constant K§ = 2K4 such that for all i, I, h and 
y G W^fT^R 3 ) one has 



Qi({h}xhSi) 



y dx 



y dx 



K 5 

< ^l|Vy|| L 2 (T h ;M3) . 



* Qi({h}xhSi) 

We now apply inequality (14. 5p to the sequence yC^) o Q x o — y c x to obtain 

||yCy) o Ch o p(M - y? - (/ y(^) o Q x o P^) dx - C Ll )\\ LH{hjM)xSi; 

J{Li}xSi 

< K 2 \\V{y^ o Q x o P(%)) - ( („?)' I I )|| L2((ft , )il)x s 1 ; R 3 ) . 



(4.8) 



Now, using the assumptions (|4.2p and (14. 3p we obtain that \\y^ o Q x o p( fe j) - yf 
The estimate (|4.7p now implies 



2 ((%,Li)xSi; 



lim ||y^ o Q x o p( fc i) - C |Ua {{ ^ }x5i;R . 

By applying (|4.8p for Z = 1 and z 7^ 1 to the sequence y\ h i) we obtain 

K 5 



0. 



(4.9) 



^ dx 



< 



||Vy(M|| 



L 2 (T h J;I 



Now we change the variables in the integrals on the left hand side (also note that T hj C £l hj ) to obtain 

K 5 



f y [h]) oQiO P^ dx - -f y^ o Q x o pfo) 
J{hj}xSi J{hj}xS! 



< 



||Vy^)|| 



L 2 (n ft J;R3)- 



Therefore (14. 1|) and (j4.9j) imply , g , y^'oQioP^) ->■ C for all i's. By applying the inequality (14.61) 



to the sequence yC 1 -*) oQ^oP^ — y| for z / 1 we obtain that HyC^oQ^oPw) — y^ ||wi.2((h-,£i)xSi;M 3 ) 
Then (|4.4p follows immediately from f|4.7|) for % = 1 and using the fact that \y\ — ColU 2 ({/i-}x5j;i 

I^I^HyC^-Coll^O. 

In the following we use the notation 



0. 



(y,d',d- i ) = ((y 1 ,di,di),...,(y n ,d 2 n ,df l )) 

to collect deformations of all rods. 

Combining the results of Theorem 13.11 and Lemma 14.11 we obtain the following result. 

Corollary 4.2 Let (yW) C W 1 ' 2 (Q h ; M. 3 ) be such that 



limsup [ dist 2 (Vy w , SO(3))dx < +00, 
h->o h Jn h 

j 

>J{Li}> 



(4.10) 



lim -f yQ 1 ^ o 

J{L!}xSi 



Qi o pfo) dx := C Ll G M 3 . 



(4.11) 



T/ien /or every sequence in R + converging to there exist a subsequence (hj) and y i G W ' ((0,Lj); 
d 2 ,d 3 G iy 1 ' 2 ((0,L i );R 3 ) suc/i that for Hi = (yj d 2 d 3 ) one has 



(y,d 2 ,d 3 ) G^l:={((y 1 ,d 2 ,d 3 ),...,(y n ,d 2 ,d 3 )) G (iy 2,2 ((0, Li); M 3 ) x iy 1,2 ((0, Li); M 3 ) x iy 1,2 ((0, Li); M 3 )) 
x ... x (I^ 2 ' 2 ((0, L n ); M 3 ) x ^' 2 ((0, L n ); M 3 ) x ^ 2 ((0, L n ); M 3 )) : 
R, G SO(3) a.e. and y x (0) = • • • = y n (0),Ri(0)Qf = • • • = Rn(0)Q^} 



and 



1 

lim -3 £ ||y^') o Q, - A(</;, d 2 , d 3 ^"~ 



(4.12) 



w/iere D i (y i ,d i ,d 3 )(xi,x 2 ,x 3 ) = y^xi) + x 2 d i (x l ) + x 3 d 3 (xi), /or x G (hj,Li) x fyS*. 



Proof. ^From (|4.10p it follows that the assumption of Theorem 13.11 is fulfilled. Therefore there exist a 
subsequence (hj) converging to and Rj G Ty 1,2 ((0, Lj), SO(3)), i = l,...,n such that Ri(0) = R 2 (0) = 
• • • = R n (0) in the sense of traces and 

1 n r — 
lim 72 E / h l|Vy (/l %) - Ri(a; • U)fdx = 0. 



We rewrite this convergence to obtain 



1 n f — 

= lim-^V/ HVy^^Qi^-R^Q^-ti)!! 2 ^ 

i n r _ 

= lim Ta E / l|V(y (,lj) o Qi)(x)Qf - R(* • e 1 )\\ 2 dx 

i n r _ 

= lim -jV / llV^oQ^-R^xOQifdx. (4.13) 



Now we define 



(y^RiCsOQiCi, y?(0) = 0, 
d? = Ri(xi)Qie 2 , 
df = Ri(xi)Qje 3 , 
Ri = Ri(xi)Qj. 

Since R G I4 7l,2 ((0, Lj), SO(3)) it follows that y° G tf/ 2 ' 2 ((0, L<); M 3 ) and R; G H /1 ' 2 ((0, Lj), SO(3)). By 
the trace property of Rj we obtain Ri(0)Qf = • • • = R n (0)Q^. 

In the sequel we want to apply the Lemma 14.11 Therefore we check the assumptions of the lemma. 
First, we estimate the norm of a matrix by the distance of the matrix to SO (3) and the norm of an arbitrary 
rotation to obtain 

f \\Vy {h ^fdx < 2 f dist 2 (Vy ( ^ ) (x),SO(3))dx + C/i 2 . 

Using (jlTOj) we obtain that flUD is fulfilled. 
Changing the coordinates in (|4.13|) we obtain 

n „ 

= lim V / ||V(y ( ^ } oQ i oP^))(x)VP ( ^ ) -R^xOQ^I 2 ^. 

This implies that (j4.2j) is fulfilled with defined above. The assumption (|4.3p is fulfilled by (|4.1ip . 
Therefore we can apply Lemma l4~T1 to obtain that for Co := Cl x — Vi(Li) we have (|4.4p and 

n 

lim Y, \\y {h3) Q< ° p(kj) ~ ViWw^UhjMxSi) = 0, (4.14) 

1=1 

where y^xi) = C + y° (a*). Since y • = (y?)' from (|4.14p and (|4.13p we obtain (|4.12|) . ^From (|4.12p and 
the estimate 

||y^) o Q, o p(^) - yi \\ L 2 {{h . }xSi) < C\\y^ o Q, o - yi \\ w ^ {{hjM)xSl ) 

(for details see the proof of Lemma I4.ip we obtain 

lim ||y^)oQ <0 P(*i)- yi \\ L 2 {{h . }xSi) = 0, i=l...,n. (4.15) 

Now, (IPj) and (j4T5l) imply 

|5i| 1/2 ||yj(/»j) - C || = - C \\L^({hj}xSi) -> 0, 

for all i = 1, . . . , n. This implies that y 1 (0) = . . . = y n (0) = Co- 
Thus we obtain that (y, d 2 , d 3 ) G A 



Remark 4.3 The structure of the functions Di(y i , df, df) defined after $4.12 ) is essentially one-dimensional. 
It stands as the limit displacement for the i-th rod. The function y i describes deformation of the middle 
curve of the i-th rod, while the vectors df and df span the normal plane of the deformed middle curve (since 
R-i = iy'i df d3i) G S0(3)j. Since the rod is assumed thin, variables x 2 and x 3 (cross-sectional coordinates 
of hSi) are of order h so the terms involving these terms can be considered as first correctors to the leading 
order approximation y i of the i-th rod. Note as well that the convergence \4-12ty will be the one which will 
be used to formulate the asymptotics of the infimizing sequence. 



Proposition 4.4 Let the functional I is defined by 
I(y,d 2 ,d 3 ) = 



n 

T- 

i=i 
+00 







qi(RjB!i)d Xl if(y,d 2 ,d 3 )eA, 
otherwise 



where Rj := (y' i: df , df) , while the class A is given by 



A:={(( yi ,d(,di),...,(y n ,di,d 6 n )) G ^((0,Ii);R i )x^((0 ) L 1 );R il )x W l >\$M) 
x ... x (W 2 > 2 ({0,L n );R 3 ) x W 1,2 ((0, L n ); IR 3 ) x V^ 2 ((0, L n ); R 3 )) : 
Ri G SO(3) a.e. and Vl (0) = ■■■ = y n (0),Ri(0)Qf = • • • = Rn(0)Qj}. 

The quadratic forms q\ : M 3 ^, — > [0, 00) are defined by 



ga(A) := min 

aeW 1 ' 2 (S i ;R 3 ) J Si 







( ^ 












X 2 




d 2 a 


d 3 aj 















where 



4(G) 



d 2 W 
dF 2 



(Qi )(GQ- ,GQ- ). 



Then the following two statements hold. 

• (liminf inequality) Let y { G VF 1,2 ((0, Li);R 3 ), df,df G L 2 ((0,L;); 
(hj) C (0,oo) converging to and every sequence (y^) C W l,2 (QJ 1 ^ 



(4.16) 



Then for every sequence 
such that 



1 n 



j^oo /l 2 . 

J 2 = 1 



0, 



(4.17) 



where Di are defined in $4.12 ), we have that 

I(y,d 2 ,d 3 )< liminf ^E^\y^); 



3 ^00 hj 



(limsup inequality) For every sequence (hj) C (0, 00) converging to and for every y i G W > ((0, Lj); 
d 2 ,d 3 G L 2 ((0,Li);lR 3 ) t/iere exist a sequence (y^)) C W 1,a (n^';R 3 ) suc/z ifarf 



1 n 

lim -2 J] ||y^ o Q, - £)j(^,d 2 ,df) 11^,2 



j->ao h 2 . 

J i=l 



{(hj,Li)xhjSi 



and 



lim ^E^\y^) = I(y,d 2 ,d 3 ). 

j->oo hj 



Remark 4.5 ^4s it is noted in Remark 3.4- in [24] each minimization problem in \4- 16\ ) has a solution 
and it can be equivalently computed on the class of functions 

Vi = {a G W h2 (Si;R 3 ) : j adx 2 dx 3 = J Vadx 2 dx 3 = 0}. 

It can be also shown that for every i the minimizer is unique in Vi and that minimizer in Vi depends linearly 
on the entries (ay) of A. Hence q\ is in fact a quadratic form of A. In the isotropic case (W(F) = W(FH) 
for every F G M 3x3 and R G SO(3)j for every i we have 53(G) = ^f(i)(G,G). In this case there are 
also some explicit formulas for q\ (see Remarks 3.5. and 3.6. in [24]) ■ 



Proof. Let us first prove the liminf inequality. 

Let (y,d 2 ,d 3 ) G A and let < hj ->■ and y( h i) C W 1 ' 2 ^; R 3 ) satisfy (HTTP . Let us also fix 5 > 0. 
Then, after rescaling each convergence in the sum (|4.17p to the fixed domain (<5, Lj) x Sj we obtain that 
for every i G {1, . . . , n} one has 

\\y^ oQ,o pCy) - y i \\ w i,2 {{6!Li)xS .. R3) -> 0, 

||(l^(y^) o Q, o P(^), ^d 3 (y^ o Q, o P^)) - (d?, df )|U 2{(5)ij)xSj;R 3 xR 3 ) -)■ 0. 

Now, by using Theorem 3.1. from [23] on each rod separately (applying it to the energy density functions 
W<% (F) := W(FQf)) we conclude that for every 5 and for every i we have 



4) 

'Li 



\ f 1 g^RfR^x < liminf-^/ ^^(/^oQ.oP^)))^ 

2 ^5 ■? ->0 ° J(S,Li)xSi 

= liminf i / ^(Vy^ ) (Q i a;)Q i Qf)^ 

j^oo ft. J(s,L i)x h n S, 



liminf i / W^Vyf^Oz;))^, 
j-Kx> ft. jQ.((s,Li)xh4Si) 



«-j JQi({S,Li)xhjSi) 

where we have used the notation = (c?i ^2 ^3). By summing all these inequalities we obtain that for 
every 5 > one has 



1 r hi 1 

£± / 4(RfR0^! < liminf -E^)(y^). 

2 J s J^OO ft^ 



By letting 5 — > we obtain 

I(y,d 2 ,d 3 ) < liminf i^^(y^)). 

j^OO ft| 

Let us now suppose that (y, d 2 , d 3 ) ^ A We have to see that for every sequence (y( hj >) C W rl ' 2 (f2 /l J';R 3 ) 

such that (|4.17p holds one has lim infj\_ +00 ^E^ h '\y^>) = +00. Let us suppose the opposite, i.e., 

3 

lim inf j-s.oo -^E^ h ^(y^ h ^) < 00. Using the property of the stored energy function W we estimate 



C * f dist 2 (Vy^)(x),SO(3))dx < * tfte)(yto>) < 00. 



(4.18) 



^From the convergence (14.170 one can easily conclude, using the continuity of the trace operator and the 
fact we can control the change of the domain (similarly as in Lemma 14.11 and Corollary 14. 2p , that 

lim f y( fc i) oQ lD p( h i) dx := y^Li). 

Thus the assumptions of Corollary 14.21 are fulfilled and we can conclude, by the uniqueness of the limit, 
that (y, d 2 , d 3 ) G A, which is a contradiction. 

To prove lim sup inequality we have to construct the appropriate sequence. Let us take (y, d 2 , d 3 ) G A. 
Let us in addition suppose y { G C 2 ([0, Li]; M 3 ), d-,d 3 G C 1 ([0, Lj]; R 3 ) (note that y { G C 2 ([0, Li]; M 3 ) is an 
immediate consequence of d 2 ,d 3 G C 1 ([0,L i ];R 3 ) and R G SO(3)). Let us define 

y( h j) in the following 

way 

y( h i\x) = yj(0) + Rj(0)Qf x, for x G (the definition is not ambiguous), 
y(M(Q,.p(M( x) ) = y , {xi _ hj) + h . y i m + ft jX2 ( d 2 (xi _ ^ _ a (^)( Xl )( d 2y (0)) 

+ft i2 ;3(d 3 (xi - fy) - a ( ^ ) (x 1 )(d 3 )'(0)) + /i 2 $(x), for x G (fy.Li) x 5,. 

where G C 1 ([ftj, Lj]; R 3 ) are such that 

a<- h *\h j ) = 0, {a^ ) )'(h j ) = l, {a^ ) )'{2h j ) = 0, 

a^^(xi) = 0, for x\ > 2hj, swp\\a^ hj ] \\oo <Chj, sup||(a^)'||oo 

< 00 



(e.g. a^)(xi) = — ^-x\ + 8x1 — 4hj, for x\ G [hj,2hj] and otherwise). The functions (3\ : 

[0,Li] x Si — > M 3 are chosen such that /3{{x) = ^(x^P^x), where 7? G C 1 ([0, L,]; M) are such that 

Tf'(xi) = 0, for X! < 7? (a*) = 1, for x x > 2hj, || 7 ^||oo < C, IKtO'IU < ?- 

(e.g. 7^'(x) = -^x 3 + Ax 2 - ^x + 5) and ft G C^fO,^] x S^R 3 ). Then we have 

3 j J 

Vy ( ^ } (x) = Ri(0)Qf, for x G T ftj ' , 
Vy^)(Q i P^)(x))Q i = Rj(xi - fc,-) - (0 |a^)(x!)(d 2 )'(0) | a^\ Xl ){d 3 )' {G))l hj < x ^ 2hj 

+h j (x 2 {d 2 i )'(x l - h j ) + x 3 (d 3 )'(x 1 - h 3 ) I d 2 (3{(x) 1 d 3 $'(x)) 
+/ l ,(-x 2 (a^))'(x 1 )K 2 )'(0) -x 3 (a^)) / (x 1 )K 3 )'(0) | | O)^.^^. 
+h 2 (d 1 p{{x) I I 0), for x G {hj,Li) x Si. 

Note that y^ G C 1 ^ ; M 3 ). It can be easily seen, by the dominated convergence theorem, that 
for every % we have linx,--^ j^-\\y [hj) o Qi - yi\\L2((h 3 ,L,)xh 3 s l ;R 3 ) = and lim j->oo j^W^U^ QiQi - 
^i|lL 2 ((/i J ,L 4 )x/i : ,s , i;iR 3 ) = which together implies that y( h ^ satisfies (|4.17p . Now we have to prove the 
limsup inequality for this sequence. 
Let us for x G (hj , Li) x Si denote 

B (M (x) = R,(xi-^-) r Vy^(Q^)(x))-Qf ^ (419) 



Then 



B) hj) (x)Qi = Ri(xi - / l ,) T (x 2 (d 2 ) / (x 1 - fy) + x 3 (d 3 )'(x! - fy) I 3 2 $(x) I 9 3 /3Kx)) 



-R i (x 1 -h j f (0 |^^iZ(d 2 )'(0) I ^^(df)'(0))U.< Xl < 2/lj 



+R t (x 1 - /i J ) T (-x 2 (a^')) / (x 1 )(d 2 )'(0) - x 3 (a^'))'(x 1 )(d 3 )'(0) | | 0)1^.< X1 < 2 ^. 
+h j -R i (x 1 -h j ) T (d l p{(x) I I 0). 

Note that for every 5 > one has 

Bf j) (x)Q l -> R 4 (x 1 ) T (x 2 (d 2 )'(x 1 ) + x 3 K 3 )'(x!) I a 2 /3 i (x) | d 3 &(x)), a.e. x G 

For every z, we look at the sequence (fj)j of functions /j : (0, Lj) xSj4 [0, +00) defined by 

/j(x) = 0, for x G (0, /ij) x Si, 

fj( x ) = -L W (Vy^)(Q 4 P^)(x))) = ^W(Qj + h j B^\x)), for x G (h^L-) x 

where the equality in the second line holds by the objectivity of W. Since W is C 2 in the neighborhood of 
SO(3) and has extreme on SO(3) and is bounded, by the Taylor theorem, for every i, one has 

fj(*) -> \qi(Kj{x 2 (d 2 )' + x 3 (d 3 )' I d 2 Pi IdaA)), a.e. x G (0, L t ) x 

Since B^ is bounded the sequence /j is also bounded in L°°((0, Li) x S^M 3 ) also by the Taylor theorem. 
Thus by the dominated convergence theorem we have 



lim i / . W(Vy^)dx = hm / \w(Vy^\QiP^\x)))dx 

= Hm f f]{x)dx = \( q l(Kj(x 2 (d 2 )' + x 3 (d 3 )' | d 2 P t |d 3 /3J). 

J ro. la x a 2 J en. l.Ix a 



Also note that for the chosen sequence y^- hj \ for every j, one has W(Vy( h ^ \ T hj (x)) = and thus 
lim i f W(Vy {h ^)dx = V lim i / , W(Vy (h ^)da; 



n 

El 



i=l 



q l(Rj(x 2 (djy + x 3 (dty\d 2 f3 l {daft)) 



(0,Li)xSi 



Thus for (y,d 2 ,d 3 ) G A s.t. y f G C 2 ([0, LJ; R 3 ), d 2 ,df G C 1 ([0, Lj]; R 3 ) and arbitrary & G C 1 ([0,.L i ] x 
Si,;~R ) we have that there exists a sequence ( y (Aj)) c W 1,a (n^';R 3 ) such that 



1 n 

lim VHyCMoQ.-A^^d: 



j->oo h 2 . 

J 1=1 



and 



lim -Ll5(Ai) („(**>) = 1 /" g » (Rf (x 2 (d 2 )' + x 3 (d 3 )' | d 2 /3, | a 3 /3J). 



Now let us now consider the general case and take an arbitrary (y, d 2 ,d 3 ) G A For every i we choose 
a sequence (R?' } ) C C 1 ([0,L i ];M 3x3 ) such that R^ -»• (y^,d 2 ,df) = R; in VF 1,2 ((0, Lj); M 3x3 ). By 
making a slight correction, namely taking Rp^ = Rj(0)R^(0) — 1 R^ ? (this can be done for j large enough 
due to Sobolev embedding theorem) we also have (R^) C C 1 ^, Li}; M 3x3 ), R^ -> Rj = (y^,d 2 ,d 3 ) 

in V^ 1 ' 2 ((0,L i );M 3x3 ) (this follows from trace theorem) and R- j) (0) = Rj(0). Now take Rj j) = ITR^ 
where II : M 3x3 — > M 3x3 is a smooth function in the neighborhood of SO (3) defining projection from the 
neighborhood of SO (3) to SO (3). We define 

f.l:i 



(J)< 



(xi) := yM + / U {j) (s)e 1 ds, d\ {3) = K^\ Xl )e k , for k = 2,3 



y 



Then((y^\dl' {j \4'^...(y^\d 2 / j \d^ {j) ))eAandy\ j) in C 2 ([0, L,];R 3 ), d 2 ' W , df W G C^O, I*];] 

and we also have that ((y\ j) )', df {j) , d 3 / j) ) = is converging to R; = (y' i: d 2 , df) in iy 1 ' 2 ((0, L 4 ); M 3x3 ). 
The functions are chosen in the following way. We choose cxi(x\,-) G V% (see Remark 14. 5p to be the 
solution of the minimum problem defining qf (Rj (xi)R.(xi)) (affine function of Hf(xi)'R' i (xi)). Now take 

/3j = RjCtj and take G C 1 ([0,Lj] x Si;R 3 ) defined by convolution (first by first variable and then by 
last two variables) such that (3^ — > and dk/3f — > dkfti (fo r k=2,3) in L 2 ($7;R 3 ). By an application of 
the Nemytsky operators theory (see [U p. 15]) we have that for every i 



g 3 ((Rp' ) ) T (x 2 (d^)'+x 3 (d 3 )' | d 2 (3Y> | d 3 PY'))dx 



j2y 



3\/ 



[0,Li] 



lO") 



qliixiRJidlY+XsRiXdi)' I d 2 a, I d 3ai ))dx. 



[0,Li] 



Therefore, we can assume (by taking a subsequence) that 



/(y,d 2 ,d 3 ^ 



1 



qi{(R { p) T (x 2 {d^>y + x 3 «' u; )' I W | 



[0,ii] 



3,0) v 



iti) 



1 

< -. 

J 



For a given j, from the previous part of the proof, we can find y^ 1 ^ G W 1 ' 2 (Q' 1 -'; R 3 ) such that 



p£||y^oQ,-A(y l w ,d 



(i) w 2 >G0 w 3 >0')mi2 



IW /1 . 2 ((/i i ,L i )x/i j Si 



1 

< - 



and 



'(0,Li)xSi 

By the triangle inequality we have that yw) satisfies (|4.17p and 



^((RfO'^^K^O' + ^K^)' I M J) I WO) 



< 



J_ E^(y^)-I(y,d 2 ,d 3 



0. 



The case (y, d 2 , d 3 ) G" .A is obvious. 



5 Minimization 



Since we can not formulate the problem of junction on a canonic domain in a simple way we have to adapt 
techniques of T-convergence and use the asymptotics of the infimizing sequence in the form (|4.17p . 

We suppose that the external body force is given by the density G L 2 (£2 h ; R 3 ) and that the external 
surface force is given by the density g { r h) G L 2 (dft h ; R 3 ) (we assume both are dead loads). As is usual in 
lower-dimensional modeling the scaling of the surface force densities is different at rod ends and the lateral 
boundary. Therefore we introduce the notation 

n (h) _ „0)i „(h) _ „(h)\ 

9 r i -9r' l8fi''\u5L 1 Qi({L i }x/i5i)) 9re ~ 9r ' lu? =1 Q i ({L j }xft5 i ) ■ 

We give the result for Neumann boundary condition on the whole domain, i.e for the pure traction problem. 
Therefore we suppose that the resultant of all forces is zero, i.e. f Qh fi h \x)dx + f dQh gi h \x)dx = 0, and 
look for the minimum that satisfies /q 1 ({l 1 } x /i5 1 ) V ( x )dx = 0. 
Theorem 5.1 For every h we define the functional 

jCO („)=/" W{Vv(x))dx- [ f ( r h \x)-v{x)dx- [ gl h) (x) ■ v(x)dx 
Jn h Jn h Jan h 

in the space 

V h = {v£ W^(Q h \ R 3 )\ / v{x)dx = 0}. 

Let the scaling of loads is as follows 

Ah) (h) (h) 

Ah) = Jr (h) = 9rj_ (h) = 9re_ 

J h 2 ' 9l h3 , 9 e h 2 » 

, in n(h) ii 1 ii (h) || j 1 II (h) n2 

where ^{P '^(o^rs), ^\\g e U^u^L^xhS^) and Ji \\g l \\ L 2( dU h\ u n =i{L . }xhS .. 
over, let us suppose that 



are bounded. More- 



[ fi h \x)dx + 
Jn h 


[ 


[x)dx 


= o, 


\f(h)( x) -f.((p(h)y 




)fdx 


= o, 


\gf\x)-g H ({P^y 




)fdx 


= o, 


gi h \x) - g ei ((pWy 




)fdx 


= o, 



(5.1) 
(5.2) 

u * — * nr tr<h 
i=i JU i 

fniEl / \\gt\x)-g u {{P^)-\QZx))\\ 2 dx = 0, (.,.:•>) 

n 

l im J^T2 \\gi h Hx)-g ei ((P^)- 1 (Qjx))\\ 2 dx = 0, (5.4) 



1=1 



'Qi({ii}xftSi) 

where f { G L 2 ((0,L;) x S f ,R 3 ), g u G L 2 ((0,L*) x dS f ,R 3 ), g ei G L 2 ({LJ x S^R 3 ), for i = 1, ■ ■ ■ ,n. 
Then we have that \ mi veV h j( h \v)\ < Ch 4 . Let us take the sequence y^ 1 ) G V h that satisfies 

J {h \y {h) )< inf J (h \v) + o(h 4 ), (5.5) 

v£V h 

(o(/i 4 ) means that lim^o ' = Let the sequence (hj) converge to 0. Then there exists a subsequence 
of (hj) (still denoted by hj) and (y,d 2 ,d 3 ) G A such that 



&pE ° Qi " D *(w < ,^,d?)||^i,. ((Ajfli4 ) xftiW) = 0. (5.6) 

3 i=l 

The limit (y,d 2 ,d 3 ) minimizes the functional 

-A r L ' r 

J(y,d 2 ,d 3 ) = I(y,d 2 ,d 3 ) - V / / f i (x)dx 2 dx 3 -y i {xi)dx l 

i=1 Jo Js, 

~X/ / / 9u{x)ds ■ yi{xi)dxi - V / g ei (Li,x 2 ,X3)dx 2 dx 3 ■ y^Li) 

~\ J0 JdSi ~1 J{Li}x-Si 



in the space V\ = {(y,d ,d ) € A: yi(Li) = 0}. Moreover, the energies converge to the energy of the limit 

lim-^jW(yW) = J(y,d\d 3 ). 

Proof. STEP 1 (a priori estimate for the total energy and y^): Let us estimate | j( h \i + a^ h ')\, where i 
is identity mapping and S IR 3 is chosen such that i + a^ 6 V h (such exists and is unique). Using 
(jlTTj) and W(T) = we obtain 

|jW(* + aW)| = | / f ( r h \x)-i(x)dx+ [ gW{x) ■ i{x)dx\ 

< Ch 3 \\f w \\ L 2(nh) + Ch 7/2 \\g[ h) \\ L 2 {dn h\^ u n =iQi{{Li}xhSt)) 

+Ch?\\gf> ||i2(uf =1 Q<({ij}xfeSi)) < C^ 4 

Then from (|5.5p we conclude that ^J^ h \y^) < C. From this we want to conclude that ^ J Qh dist 2 (V?/^, SO(3)) 2 (fx < 
oo, so we have to estimate the energy from below 

^jW(yW) > C w ^^dist 2 (VyW,SO(3)) 

1 1 

-^2 II/ w IIl 2 (^)II?/ (/i) IIl 2 (^) - -^\\9i llL2( a Qh\ (u ^ i{L j xftSi) ||y (/i) || L 2 (9n h\ (U n =i Q i({L j x/i5i)) 

-^2llfl , e h) ||L 2 (U^ 1 Q l ({i 1 }x/ l 5 l ))ll?/ (?t) llL 2 (U^ 1 Q l ({^}x/i5 l )) (5-7) 

C w ^^dist 2 (VyW,SO(3))-C(i||^)|| L2(nh) + ^||yWl 

+ ^l|y W ||L 2 (U™ =1 Q,({LJx/ i S l )))- 



> 



Mdf^VU^Q^LJx/iSi)) 



We have 



\y {h) \\ L im = J2 \\y {h) hi(c?) + \\y {h) \\Li(T») 



i=l 



In the same way as in Lemma [4.1l we conclude that there exists constant C independent of h (using rescaling 
aS^'fai, X2, X3) = {x\,hx2,hx3)) such that for every i and every y 6 W 1,2 {Q h ; IR 3 ) we have 



y — X ydx 

Qi({Li}xhSi) 



y — X 

Qi({h}xhSi) 



L 2 (C* 



< C||\7 2/11 L^C^M 3 )' 



< C\\Vy\\ L2(ch 



(C'';R3)- 



(5.8) 
(5.9) 



^From this we conclude that there exists a constant C independent of h such that 

C 



ydx — -r ydx 

Qi({h}xhSi) J QidLiyxhSi) 



< tHIv</II L 2 (c ,. 



(5.10) 



By using scaling aS h '(xi, X2, x%) = (hxi,hx2,hx^) we conclude that there exists C independent of h such 
that for every i, I 



y 



y dx 



Qi({h}xhSi) 



< hC\\Vy\\ L2 



L 2 (T h ;R 3 ) 



y dx 



Qi({h}xhSi) 



y dx 



QMhjxhSi) 



< 



c 



(5.11) 
(5.12) 



Using estimate (J5.8P for i = 1 and the fact that € V h we conclude that 



y {h) \\mc^ < C||Vy w || L2 ^ :K3) < C(||dist(VyW,SO(3))|| i2 ^ ;K3 , + h). 



.00 



(h) 



Using estimate f|5. 10|) we conclude that 

II / y {h) \\ < j\\W h) \\ L >(c^) < c(h <ast(W h \so(3))\\L*(c?m + 1)- 

Using estimate (|5.11[) we conclude 

\\y {h) \\m^m < ^ 3/2 ll / y {h) II + ch\\v y ^\\ LHTh . R s y 

Since 

W^h^m < C(|| dist(VyW, SO(3)|| L2(nh . R3) + /i 3/2 ), (5.13) 

we conclude 

llv (ft) [|r»cz*;R3) < C(/ l 1 / 2 ||di S t(VyW,SO(3)|| L2(nh;M 3 ) + h 3 / 2 ). (5.14) 
Using estimate (j5.9j) and ()5. 12[) for I = 1 we conclude that for every i 

\\y {h) \\ LHc ^ m < C(||dist(VyW,SO(3))|| L2(nh;R 3 ) + h). 

Thus we have 

\\y {h) \\mn^) < C(||dist(VyW,SO(3))|| L2(nh;R 3 ) + h). (5.15) 

In the same way one can analyze traces. First we start from the trace inequality on the cylinder Q = 
(0, 1) x Si. For every y £ W 1,2 (Cj;R 3 ) we have that there exists constant C such that 

\\y-f yWmdCA < C\\y - f vllwi.arcfcRS) < C||Vy|| L2 (c i; M3), (5.16) 

J{l}xS, J{l}xS t 

f y\\L2(dCi) < C\\y - f y\\w^(cvm ^ C|I v 2/IIl2(c 8 ;]R3)- (5.17) 



12/ 



By using appropriate scaling and rotation we have that there exists constant C such that for every and 
y G PU 1,2 (C^;M 3 ) we have 



112/ — T ylU^QittMx/^);]^) < C 'l|Vl/|| i 2 (C h. R 3 ) , (5.18) 

jQi({h}xhSi) 

\y~f l/llL2(Q i ({L^x^5 i );R3) < C|| Vy || L 2( C h ;R 3) , (5.19) 

■/Qi({ij}xA5i) 

111/ — T l/ll^fQirwxhSOjRS) < C , ||Vy|| £a(0 rfc. H s ) , (5.20) 

J Qi({h}xhSi) 

\\y - T l/llL^fQiQMxhSi)^) < C||Vy|| L 2( C h ;M 3) (5.21) 



and 



I!/ - / y|U 2 (Q i ((ft,i i )xa/ l S i );R3) < CTI72 H V J/IIl 2 (C^;M 3 )' ( 5 - 22 ) 



12/ -T y\\LZ(Q i ((h,Li)xdhS i );W>) < ^71/2 II V 2/IL 2 (C^;R3): ( 5 - 23 ) 

yQ i ({L i }xfeS<) « ' 

II2/-T J/|lL2(Q i ((/i,L i )x8/iS i );R3) < C TTB H V ^Hi 2 (ChR3)) ( 5 - 24 ) 



1 

2/||L 2 fQ,(f/i.L,)xd/iS,);IR 3 ) < ^7 

'Qi({Li}xftSi) 

In the same way we conclude 



y- 7^ /fr i ^ o ^V\\l?{Cl i {{h,Li)xdhS i )-SP) < CT-iT^IIVyll^^h.KS). (5.25) 



\y — T y II 2,2(9^*3) < C/i 1/2 ||Vy|| L2{T h R3) (5.26) 



Now by using y^ G V h we have from (|5,19p and (|5.23p that 

\\v {h Hi?W{L$xhSr)&) < <?l|VyW|| x2(chK3) <C(||dist(VyW,SO(3))|| L 2 (nh . R 3 ) +/ l ), (5.27) 
H?/ (/l) |lL 2 (Qi((/i,Li)xahSi);R3) < C '^I72ll V 2/ ( ' l) llL 2 (Cf;M3 ) < C(j^\\dist(Vy( h \SO(3))\\ L2{nh .^ ) +0.^) 
^From (|5.18|) and (|5.19p we conclude 

y^<j\\W h) \\ LH ct^y ( 5 - 29 ) 



iQ!{{h}xhSi) 
^From this and ()5. 12j) we conclude for every i 

II / y {k) \\ < j\\Vy ih) \\Li iC t,RZ) + ^\\VV W h*<Tk-R*) < f IIVy^lLa^). (5.30) 

^From (15.221) and (j5.15[) we conclude 

Wy^hHQ^h^xdhs^m < c , (t T WiIVi/ (/i) ||l 2 (^;r3 ) + / 1 1 / 2 ||/ y ^\\) 

( lVy^|U 2( ^ ;R3) < C(-^||dist(VyW J SO(3))|| £8(nh;R 3 ) + /iYS)31) 



< 



^1/2 II"" IW'(U»;K*> ^~K hl/2 

^From Qgjggp for i = 1, ([Qg]) and ([533]) we conclude 

llv (ft) lli*(erW<^^ (5.32) 
^From (|5.18p and (|5.30[) we conclude for every i 

<C(||dist(V»W J SO(3))|| ia(nhiE a ) +/»). (5.33) 

^From (IQTJ) . (I5T32|) and (Q3l) we conclude 

ll^llL^a^Vu^Q.H^Ix^)) < C(^||dist(VyW,SO(3))|| i2( ^ ;M3) + / l 1/2 ), (5.34) 
Hw (fc) llL»(u r=1 Q 4 ({L i }xhS i )) < C(||dist(VyW,SO(3))|| i2(n , l;R 3 ) +/ l ). (5.35) 
By using (|5.7p , (|5.15p , (|5.34|> and (I5.35|) we conclude that there exists C 2 , C3 such that 

C w ± J^dist 2 (Vy^(x),SO(3))dx-C 2 (±- ^ dist 2 (VyW(a;), SO(3))cfc) 1/2 + 1^ < ^jW(yW) < C 3 

(5.36) 

Using the fact that for /i < 1 

dist 2 (VyW(x),SO(3))^ 7 < ^ ^ dist 2 (Vy W (x),SO(3))(ix X 7 
we conclude from (|5.36p that 

CV« 2 — C*2a < C3. 
which implies that a 2 is bounded, i.e, there exists C > such that 



i / dist 2 (Vy^(x),SO(3))dx < C (5.37) 
fr 4 Jnh 



which implies that the left hand side of (|5.36p is bounded as well. This implies 

I inf J( h \v)\ < Ch*. 
vev h 



STEP 2 (the convergence proof for and the scaled total energy): The estimate (|5.37p implies 
that the assumptions of Theorem 13.11 (compactness theorem) are fulfilled. Therefore the assumptions of 
Corollary 14.21 are fulfilled as well (with Cj_, x =0). Therefore we conclude that for every sequence hj there 
exists a subsequence (still denoted by hj) and (y, d 2 , d 3 ) = ((j/i, df, df ), . . . , (y n , d^, d 3 )) 6 A and 



n 



L m ^ E ° Q< - A(^,^ 2 ,df)|| 2 wl , 2((?t . iLi)></lj5i;M;!) = 0. (5.38) 

J 1=1 

^,From this convergence it is obvious that y l (Li) = ]imj-+ 00 /q 1 ({£, 1 } X /i-Si) H^ hj \x)dx = 0. Thus we have 
proved that (y,d 2 ,d 3 ) G V/ and what is left to prove is that it minimizes the functional J in V\. We can 
use the standard argument from T-convergence, although we have variable domains (and can not apply 
the T-convergence directly). Let (y a ,d 2 ,d 3 ) G V\ and (y a ,d 2 ,d 3 ) / (y,d 2 ,d 3 ). We have to prove that 
J(y, d 2 , d 3 ) < J(y a , d 2 , d 3 ). ^From the liminf inequality from Proposition 14.41 we conclude 

I(y,d 2 ,d 3 ) < liminf ^E^\y^). (5.39) 
By using ([51i|l . ([Q2]) and ([QT]) we have 

da; (5.41) 



"4 / „ g^^-y^^) 

*i JdT h i\LI?_AhAxhjSi 



<— Wa^W h Wv^W h < C(h {5/2) +h 2 )^ n 

hj" 1 11 l 2 (9t fe j\utL 1 {h J -}xhiS , i)ll* / Ili^arS-yu^j^jxftjSj) - ^v"? -i-^-t-u. 

^From this and (^-([53]), doTTBI) . (Q4"|) . (i5\35l) . (i5\38l and (i5\39l) we conclude that 



6 



J{y,d 2 ,d 3 ) < liminf -L jCy)( y ( h i)). (5.42) 



Let us by Zimsup-inequality from Proposition 14.41 choose t/i^ such that 



1 n 

fe, 7^2 E ° Q< " Afoi,*' d i,a> <a)ll H /i, 2 ((/ lj ,L l )x^S l ;R3) = (5.43) 



j-5-oo 



J 8=1 

and 

lim ±E^\yi h > ) ) = I(y a ,d 2 a ,dl). (5.44) 

3 

Let us choose c^') £ R 3 such that = y„ + G . From the convergence (I5.43j) we conclude 
that Um^oo c^') = 0. Thus we have that (15. 43ft is also satisfied for the sequence . We also see that 
()5.44p is also satisfied for z^ j \ Therefore, using the lower bound on W, it follows that there exists a 
constant C such that sup^ ^ J Qhj dist 2 (Vzi hj) (x),SO(3))dx < C. In the same way as before we conclude 



J(y a ,dl,dl)= lim -Lj(M(^M). (5.45) 



Finally from (|5"3j) . (|5.42p and (I5.45j) we have 

J(y,d 2 ,d 3 ) < liminf -LjO^M) < liminf 1 J^\y^ j) ) = J(y a ,d 2 a ,dl) 

j— too n*; j— >-oo tij 



3 J 3 



That the energies converge can be easily seen by standard argument in T-convergence (we first take sequence 
such that Aj( h i)(Z( h i)) J(y,d 2 ,d 3 ) and then by using conclude that lim,-^ Xj^Q,^) = 

h 3 3 

linij^oo j^J^ hi \y^ hj ^)- Since this can be done for arbitrary sequence, we have the claim). 



Remark 5.2 Using J5.^0| J and \5.J± 1\ ) by straightforward calculation from \5. 1}) we obtain 

Y ( / ( / fi{x)dx 2 dx 3 + / g u (x)ds] + / g ei {U, x 2 , x 3 )dx 2 dx 3 ) =0. 
\Jo \JSi JdSi J J Si J 

This means that in the limit model the total force is zero as well. 

Remark 5.3 Adding a constant to the solution of a pure traction problem gives a solution again, i.e., 
the set of solutions is closed under translations. Therefore, we had to control behavior of this constant in 
three-dimensional problem in order to obtain the limit. We did it by requesting that the mean deformation 
at the end of the first rod (indexed by 1) vanishes. As expected, this constraint results in the limit model 
in the constraint that the end of the first rod is fixed in the origin (y^Li) = 0). In the limit model we can 
also consider this constraint as the one which just fixes the translation since again the set of solutions of 
the pure traction problem is closed under translations. 



6 Differential formulation of the model 

In this section we formulate the weak and the differential formulation of the model. It enables us to give the 
interpretation of the limit model as the model of one-dimensional rods with the transmission conditions 
at the junction point, see (|6.6j) — (|6. 12|) . 
Let us define 



fi(xi) = / fi(x)dx 2 dx 3 + gu(x)ds, 
J St JdSi 

Fi = g ei (Li,x 2 ,x 3 )dx 2 dx 3 . 

J Si 



'Si 

Then the total energy of the limit model is given by 

" / 1 r L * r L ' - \ 

J(y, d 2 , d 3 ) = J2 [2 J Q 4(R-fR-D^i - J o • Vidxi - Fi ■ yiU)) , 

for (y,d 2 ,d 3 ) G A and +00 otherwise. 

First, performing partial integration in the force terms in the total energy functional, similarly as in 
[25], we remove appearance of y i from the energy functional. In order to do that let us denote 

Pi(xi) = fi(z)dz + Fi, i = l,...,n (6.1) 

J X\ 

and note that the force equilibrium, according to Remark 15.21 can be expressed by 

n 

J2pM = 0. (6.2) 



Then 



i=i 

n „£. n n n 

J2 Pi-Vidxi = -ViiLi) - Vpi(0) -yj(0) - V / Pi-ytdx! 

i=i J * *=i i=i j=i J ° 

n n n „^ 

= ■ - Ea(°) • yM + Y £ • y^ 

i=l i=l i=l •'° 

n n . L . 

= Y Fi ' y^ L ^ + Y J y i dxi > 



i=l i=l 



since in A deformations satisfy y\(0) = •••y n (0). Thus the total energy functional can be expressed in 
terms of Rj, i = 1, . . . , n only, by 

n fl fLi rLi 

J(R) := J(y, d 2 , d 3 ) = Y[2j Q <i(^I^)dxi - Pi- Rieid Xl 



where we have used the notation R = (Ri, . . . R n )- Thus we can split the minimization of J in two steps. 
In the first step we minimize J in the space 

TZ := {R = (Ri, . . . ,TLn) £ W x ' 2 ((0, L x ); SO(3)) x . . . x W x ' 2 ((0, L n ); SO(3)) : 
Ri(0)Qf = --- = Rn(0)Q^} 

and in the second step we determine deformations i = 1, . . . , n from the equations 

y'i = R, i e 1 , yi (0)=y , i = l,...,n (6.3) 

where the constant vector y £ R is freely determined from an additional constraint (e.g., y 1 (Li) = 0, see 
Remark 15 .3 j) . 

Now we want to find the weak formulation of the problem 

min J(R). 

Reft 

First note that Rf R^ are a.e. anti-symmetric matrices. Therefore they possess axial vectors = Sj(Rj) £ 
L 2 ((0,Li);K 3 ), i.e., ' 

R-fR-i = A Sl , i = l,...,n, 

where the notation A s stands for the matrix such that A s x = s x x. Since q\ are quadratic forms of 
the elements of Rf R^ there are positive definite matrices Hj (positive definiteness of the matrices 
follows from the fact that the second derivative of W is greater or equal to and equal to exactly on 
antisymmetric matrices) such that 

4(RfR<) = n iSl • ai . 

Thus the total energy functional can be written by 

J( R ) = Yl{\J H * S *( R • s i (K i )dx 1 - J pi ■ TLie^x^J . 

In order to obtain the weak and differential formulation of the model we need to find the Gatoux derivative 
of the functional J over TZ. Let R £ 1Z, e > and V{ £ C°°([0, Li]; M 3 ), i = 1, . . . ,n. Let us choose a 
perturbation R E = (Rf , . . . , Rfj) £ TZ of R in the following form 

Rf = e eA *'iRj ! i = l,...,n. 

In order R e to be in TZ one only needs to fulfill that 

Rf(0)Qf = e eA ^>R i (0)Qf 

is independent of i. Since Rj(0)Qf is independent of i, Rf (0)Qf is independent of i if and only if 

»i(0) = ■■• = «„(0). 

Thus in the sequel we take 

v£K t = {v = (v 1 ,---,v n ) eC°°([0, L{]; R 3 ) x • • • x C°°([0, L n }; R 3 ) : Vl (0) = • • • = « n (0)}. 
Next we need to compute the axial vectors sf of (Rf) T (Rf )'. 

(Rf ) T (Rf )' = Kfe- £A ^(e A <Hi + e £A ^B!i) 

= Rf (I - eA v . + 0(e 2 ))(eA < + 0(e 2 ))R; + RfR- 

= Rf R' + eRf A^/Rj + 0(e 2 ). 
Since Rf A,/RjX = Rfuj x x (as Rj(xi) G SO(3) a.e.) we obtain 

s £ = Si + eKjv'i + 0{e 2 ), i = l,...,n. 



Now, we plug this perturbation into the functional J: 



J(Rf ) =it,(^f ' U ^ Si + m ? v 'i) ■ ^ + sKfv'^dx! - jT ' ^ • (I + eA Vi )YL iei dx^j + 0(e 2 



Thus the stationary point of the functional J satisfies 

"Li f'Li 



n 

£ 

i=i 



/Li rLi \ 

HiSi ■ Rjv'idxi - J A Vi Rieidxij = 0, v £ TZ t 



i.e. 



itji^f R^iSi " v 'idxx - j Vi ■ R^ei x Pidx^j = 0, veTZ t . 
Thus by partial integration on every rod we obtain the equation and the boundary condition 

(RiHisO 7 + R iei x ^ = 0, Rj(-Li)HjSj(-Li) = 0. (6.4) 
Moreover, since v £ TZt, we obtain just one condition in the junction point 

n 

^Ri(0)H iSi (0) = 0. (6.5) 
i=X 

Let us now denote 

Si — R^2? Qi — R^H^Ri s^. 



Then the problem given by (|6.ip . (|6.3p . (|6.4p . (|6.2p and (|6.5p can be formulated by 

# + / 4 = 0, p i (L i ) = F i , i = l,...,n, (6.6) 

+ Rid x ^ = 0, q i (L i ) = 0, i = l,...,n, (6.7) 

Qi = RiHjRf i = l,...,n, (6.8) 

Ri = Ag.Rj, i = l,...,n, (6.9) 

t/- = Rjei, z = l,...,n, (6.10) 



^p,(0) = 0, ^(0) = 0, (6.11) 

i=l i=l 

Rl( )QT = ■ • • = R„(0)Q^, yi (0) = ... = y n (0). (6.12) 



The first five equations (|6,6p - (|6.10p are the equilibrium equations of the nonlinear inextensible rod model, 
see [24] for the derivation of the model from the threedimensional nonlinear elasticity and [4] for the direct 
foundation of the theory of nonlinear rods. See also [T7] or [30] for the rod model obtained by linearization 
of the present one. The model is written as a first order system of ODE's. Introduced unknowns Pi 
and (Ji are the contact force and contact couple corresponding to the i-th rod. The equations (|6.6p and 
(|6.7p are equilibrium equations together with the boundary conditions. (|6,8p is the constitutional law, 
(|6.9p and (I6.10p are material restrictions of unshearability and inextensibility. The conditions (|6.1ip - (l6,12l) 
are conditions at the junction. The two conditions in (|6.1ip are the equilibrium conditions and say that 
the sum of all contact forces and couples in the junction are 0. The conditions in (I6.12|) are continuity 
conditions. The first one say that the rotation of the cross-section in the junction is the same looking from 
all rods. Note here the difference between Rj and R^Q^. The matrix Rj(0) gives actual position of the 
tangent vector Rj(0)ei and the cross-section (spanned by Rj(0)e2, Rj(0)e3). Rj(0)Qjf is the rotation of 
the cross-section "in the junction" of the rod for the i-th rod (the "difference" between undeformed Qi and 
deformed R«(0) configuration). The second equation in (|6.12p say that the deformation in the junction 
point is the same for all rods. 

Thus we conclude that junction (transmission) conditions for junction of rods are given by the equi- 
librium of contact forces and couples as well as by continuity of the deformations and rotation of the 
junction. 



Remark 6.1 The minimization problem for the total energy J on A from Theorem I5.il has at least one 
solution. Thus y { G W 2 ' 2 ((0, Li); M 3 ), R; G IF 1,2 ((0, Li); S0(3)). From the differential formulation for each 
rod we can conclude a certain regularity result. For fi G L 2 ((0, Lj); M 3 ) one has that pi G VF 1,2 ((0, Z^); R 3 ). 
Therefore R i e 1 x p l G VF 1,2 ((0, Lj); ]R 3 ) as we//, so ^ G IF 2 ' 2 ((0, Li); M 3 ). f/smg i TO)) to implies S, G 
VF 1 ' 2 ((0,L i );M 3 ). Now using < fO|) we ootoi to R; G VF 2 ' 2 ((0, Lj); SO(3)). Now going back to < fO|) 
we obtain that Si G VF 2 ' 2 ((0, Lj); M 3 ), w/uc/i again using \6.9\) implies that Rj G W 3 ' 2 ((0, Li); SO(3)) and 
Vi G VF 4 ' 2 ((0,L;);IR 3 ). T/iis is the most that can be concluded for L 2 loads in this fashion. 
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